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1. Introduction 

The charm quark contributions to the nucleon structure functions are large in the small x region. 
To obtain a correct theoretical description of the structure function F2(x,Q 2 ), which is measured 
at an accuracy of 1% in this region, 3-loop (NNLO) corrections to the heavy flavor Wilson coeffi- 
cients are needed. Present extractions of the strong coupling constant from deep inelastic data [1] 
and precision measurements of parton densities require this accuracy, with important consequences 
for precision measurements at Tevatron and LHC [2]. In previous calculations the Wilson coeffi- 
cients were obtained to 0(a 2 ) [3] l . If one considers the region Q 2 /m 2 >\0 the massive Wil- 
son coefficients can be obtained through convolutions [5,6] of massive operator matrix elements 
(OMEs) [7-13] and the massless Wilson coefficients [14]. Other massive OMEs play a role in 
higher order QED correction, [15]. Based on the 2-loop massive operator matrix elements, which 
are known for general values of N up to the linear term in the dimensional parameter e [8, 10, 11] 
and the massless 3-loop Wilson coefficients [14] one may compute all logarithmic contributions 
at 3-loop order [16]. In the kinematic region of HERA, however, the logarithmic terms are not 
dominant over the yet unknown constant terms for general values of N, as demonstrated for the 
moments calculated in [6]. The O(a^) corrections to the structure function Fl(x, Q 2 ) in the asymp- 
totic region were computed in [6, 17]. Heavy flavor corrections for charged current reactions in 
Mellin space were recently given in [18], correcting some errors in the literature. 

The heavy flavor contributions to the structure function F2(x,Q 2 ) at O(af) in case of one 
massive quark are described by five massive Wilson coefficients ,L p q s ,L s g ,H? s ,H* [6] in the 
asymptotic region. Two of these Wilson coefficients, L^ s and haven been computed com- 
pletely for general values of N in [12], cf. also [16]. In [12] the contributions of the color fac- 
tors 0{T 2 N f C A)F ) to the Wilson coefficients L^ S ,H PS ,H^ were also calculated. The corresponding 
Feynman diagrams consist of graphs with a massless and a massive internal fermion line. After 
applying algebraic relations [19] these contributions to the massive Wilson coefficients can be rep- 
resented in terms of the known w = 4 set of harmonic sums [20]. The calculations of the massive 
OMEs at general value of N were performed using tf orm [21] and C. Schneider's summation 
package Sigma [22], which has received continuous updates during the last years. 

In this note we report on further progress of the 3-loop calculation of massive OMEs on contri- 
butions due to charm and bottom lines in single diagrams, massive gluonic OMEs, massive 3-loop 
ladder graphs, and cyclotomic harmonic sums and polylogarithms, which emerge in intermediate 
steps of present calculations. 

2. Two massive quarks of unequal mass 

Beginning at 0(q s 3 ), graphs with internal massive fermion lines carrying unequal finite masses 
contribute. Due to the mass ratio x = m^/m^ ~ 1/10 for the case of charm and bottom quarks, one 
may expand the corresponding diagrams using this parameter. The results for the moments N = 2 
and N = 4 of the gluonic operator matrix element AQ g were given in [13] before, extending the 
code qexp [23] to higher moments applying projectors similar to those used in [6]. With rising 

'Precise numerical implementations in Mellin space were given in [4]. 
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values of N the calculations become more and more time consuming even using tf orm [21]. Here 
we present the moment N = 6 which requested several months of run time : 
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(2.1) 



These contributions to the massive Wilson coefficients can be uniquely calculated in the fixed 
flavor number scheme with three massless quarks in the initial state. However they cannot be 
attributed either to the charm- or bottom distribution in a variable flavor scheme. 2 Despite all these 
contributions are universal, the so-called variable flavor scheme ends at the 2-loop level. In general 
we remark that the representation of the heavy flavor Wilson coefficients to O(a^) are given in the 

2 It was shown in [24] that the matching scale is rather process dependent and may differ significantly from the 
heavy quark mass mg, as comparisons to complete calculations show. This aspect, however, is often ignored. 
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fixed flavor number scheme. As has been shown e.g. in Refs. [25,26] this choice is sufficient for 
the kinematic range at HERA. 



3. Contributions due to 3-Loop Ladder Graphs 

A next topology class consists of ladder graphs containing up to six massive lines and local operator 
insertions. The calculation of these graphs is underway [27]. Here we present, as an example, the 
result for the scalar Diagram shown in the following Figure : 
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4. Contributions to Gluonic Matrix Elements 

First results for general values of N have been obtained for the gluonic OMEs of 0(NfTpCA,F- 

They are of importance to define the variable flavor scheme for the single heavy mass contri- 
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butions, cf. [6, 11,28]. As an example we show the renormalized result for A k ^q{N), cf. [29]: 
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5. Cyclotomic Harmonic Sums and Harmonic Polylogarithms 

In various calculations of the massive OMEs harmonic sums emerged, which are associated to cy- 
clotomic harmonic polylogarithms. It appeared therefore as useful to generalize the usual harmonic 
polylogarithms [30] based on the alphabet {1/jc, 1/(1 — x), 1/(1 +x)} to that of the alphabet 

*€N+,0 </<?(*)), (5.1) 
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with the kth cyclotomic polynomial and <p Euler's totient function, cf. [31]. The cyclotomic 
harmonic polynomials are obtained as iterated integrals of the letters in 21. The corresponding 
harmonic sums 3 iterate the letters of the type 

(-1)* 

n>l,m <«,/>!, (5.2) 



(nk + m) 1 

which can also be related to usual harmonic sums at numerator weights being nth roots of unity. 
The properties of these harmonic polylogarithms, harmonic sums and associated generalizations of 
the multiple zeta values [33], as well as the relations of these quantities, are given in [31]. They are 
implemented in the package HarmonicSums [34]. 
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